Robustness of topologically protected surface states in layering of Bi2Te3 thin films 
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Bulk Bi2Te3 is known to be a topological insulator. We investigate surface states of Bi2Te3(lll) 
thin films using density-functional theory including spin-orbit coupling. We construct a method 
to unambiguously identify surface states of thin film topological insulators. Applying this method 
for one to six quintuple layers of Bi2Te3, we find that the topological nature of the surface states 
remains robust with the film thickness and that the films of three or more quintuple layers have 
topologically non-trivial or protected surface states, in agreement with recent experiments. 

PACS numbers: 73.20.At, 73.20.-r, 71.20.Nr, 71.15.Mb 
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Recently, topological insulators (TIs) with time- 
reversal symmetry have attracted attention due to their 
topologically protected states \^]^. In three dimensions, 
TIs differ from band insulators in that a bulk energy 
gap opens due to strong spin-orbit coupling (SOC) with 
metallic surface states in the bulk energy gap. Several 
bulk bismuth-based alloys were discovered to be three- 
dimensional TIs '2h13|. For ordinary metals, surface 
states are so fragile that they can be destroyed by a small 
number of impurities or disorder at the surface. How- 
ever, the surface states of TIs are topologically protected 
in that impurities preserving time-reversal symmetry can 
neither destroy nor impact the topological nature of the 
surface states. 

Three-dimensional TIs are classified according to a 
topological invariant, the Z2 invariant vq [1|. Strong 
(weak) TIs have vn = 1 {vq = 0). Recent first-principles 
calculations [a, [iJ] show that bulk Bi2Te3, Bi2Se3, and 
Sb2Te3 alloys are strong TIs with a single Dirac cone 
below the Fermi level, i?/, at F (fc = 0). This feature 
was confirmed by angle-resolved photoemission spectra 
(ARPES) experiments [i,0-[n|- 

Thin films offer valuable probes of TIs as well as po- 
tential device applications. It is, therefore, important to 
identify surface states and their topological properties. 
For example, thin TI films were proposed to be efficient 
thermoelectric devices that exploit the interaction be- 
tween top and bottom surface states ,1^]. Thin films also 
have considerable advantages in a direct measurement of 
transport properties of the surface states by emphasizing 
surface states over bulk states. Very recently, Bi2Te3 and 
Bi2Se3 thin films with a thickness of a few nm were exper- 
imentally realized 16l42C|. and topological properties of 
such films have been examined [18-23]. One theory sug- 
gests that the quantum spin Hall phase of a thin TI film 
oscillates between topologically trivial {i/q = 0) and non- 
trivial (I'o = 1) states with the film thickness [2l|, while 
another study implies a topological quantum phase tran- 
sition with an oscillation in an energy gap A with the 
thickness [23[. However, ARPES measurements 18l420l| 
on thin films did not show any oscillation in either i/q 
or A with thickness. These discrepancies cast doubt on 



the robustness of topological surface states and how to 
identify them. 

In this work, we construct a method to unambiguously 
identify topologically protected surface states for thin TI 
films. Our method based on density-functional theory 
(DFT) with SOC, can be used to determine the topolog- 
ical nature of the surface states in thin TI films. As an 
example, we apply this method to Bi2Te3 thin films at 
six different thicknesses. We show that the topological 
nature of the surface states remains robust with the film 
thickness and that the surface states are topologically 
non-trivial {vq = 1) for films three or more quintuple lay- 
ers (QLs) thick. This differs from previous calculations 
[2l|, yet agrees with experiment [l9|. 

We begin with a review of the bulk properties of 



Bi2Te3 [24]. A unit cell of bulk Bi2Te3 consists of 



a rhombohedral structure with five inequivalent atoms: 
two Bi and three Te. For the lattice constants 
(a=4.386 A and c=30.497 A) and the positions of the 
five atoms, we use experimental data [25|. We calcu- 
late the electronic structure of bulk Bi2Te3 using a DFT 
code, VASP [261, within the Perdew-Burke-Ernzerhopf 
(PBE) generalized-gradient approximation (GGA) [27| . 
Projector-augmented-wave (PAW) pseudopotentials are 
used [28(. Plane waves with a kinetic energy cutoff Ec 
of 175 eV are used as basis sets and 146 irreducible k- 
points (iVfe = 146) are sampled. Figure [1] shows the bulk 
band structure computed without and with SOC at time- 
reversal invariant momenta, F, Z, F, and L, with Ef = 0. 
At the Dirac point (F), SOC inverts the conduction and 
valence bands with opposite parities [5(. 

Let us discuss the structure of a Bi2Te3(lll) film and 
parameter values for band-structure calculations of the 
films. Along the (111) direction (trigonal axis), a Bi2Te3 
slab is built in units of a QL, which consists of two Bi and 
three Te layers that alternate [inset in Fig. EJb)] . In each 
atomic layer, the Bi or Te atoms form a triangular lat- 
tice, and their in-plane positions coincide with those in a 
(111) surface of a faced-centered-cubic lattice. For each 
QL [inset in Fig. [2i;b)], the Te(l) and Bi(l)/Bi(2) lay- 
ers bond more strongly than the Bi(l)/Bi(2) and Te(2) 
layers, since zi < Z2- Neighboring QLs are separated by 




FIG. 1: Band structure of bulk Bi2Te3 (a) without SOC and 
(b) with SOC at symmetry points F, Z, F, and L. 



TABLE I: The number of crossings A'^c of the surface states 
across Ef, the energy gap A at F, and the indirect energy gap 
A-Eind (eV) as a function of slab thickness (nm). Ec (eV) and 
Nk are the energy cutoff and the number of k points sampled, 
respectively. 



no of QLs 


thickness 


Ec 


Nk 


A (eV) 


N, 


ABind 


1 


1.0166 


300 


48 


0.4338 





0.301 


2 


2.0332 


300 


96 


0.1319 





0.057 


3 


3.0497 


500 


341 


0.0261 


1 


No gap 


4 


4.0663 


500 


341 


0.0070 


1 


No gap 


5 


5.0829 


500 


341 


0.0090 


1 


No gap 


6 


6.0993 


500 


341 


0.0055 


1 


No gap 



^3 (> Z2), which causes the Te(l) layers between neigh- 
boring QLs to interact via weak van der Waals forces, 
allowing exfoliation 11 71 . For the slab electronic struc- 
ture, we use VASP [2g with the PBE GGA and PAW 
pseudopotentials as in the bulk case. We do not per- 
form reconstruction of the surfaces. Compared to the 
bulk case, the electronic structure of the slabs, especially 
of surface bands, converges more slowly with parameter 
values {Ec, N^, and the number of vacuum layers). For 
the six slabs with different thicknesses considered, a large 
vacuum layer equivalent to 5 QLs (=50.829 A) is added. 
The values of Ec and Nk used for the six slabs are listed 
in Table m 

Self-consistent DFT calculations with SOC are first 
carried out using the parameter values listed in Table |T] 
until the total energy converges to within 1 x 10^^ eV. 
The two-dimensional band structure is then computed 
non-self-consistently using a charge density distribution 
obtained from the previous self-consistent calculations. 
Accuracy in the structure of the surface bands depends 
on the accuracy in the charge density distribution. Thus, 
extremely well-converged self-consistent calculations are 
required for an accurate identification of the surface 
bands. 

We now introduce our method to identify surface states 
from the slab band structure. Surface bands refer to 
bands localized on the top or bottom surface layers of a 
slab [29,], but not all bands in a slab belong to the sur- 



face bands. To identify surface bands, a wave function 
at a given energy band and fcii (momentum parallel to 
the surface) is projected onto spherical harmonics which 
are non-zero within some radius around each ion. Then 
two criteria are applied to the wave function projections. 
In criterion 1, the surface bands are identified based on 
a critical percentage of the projections onto the top two 
or the bottom two atomic layers. The critical percent- 
age depends on the slab thickness. In criterion 2, the 
surface bands are identified according to a critical per- 
centage of the projections onto the top or the bottom 
QL. Criterion 2 is suitable for multiple QLs. The den- 
sity of states (DOS) projected onto each atomic layer of 
the topmost and bottommost QL reveals that the DOS 
projected onto the inner Te(l) layers more closely resem- 
bles that projected onto the outer Te(l) layers than those 
onto the middle Te(2) layers [inset in Fig. [2l[b)]. Addi- 
tionally, the electron density spreads over all five atomic 
layers within the topmost or bottommost QL. 

Band structures of Bi2Te3(lll) slabs of 1-6 QLs, are 
shown in Fig. [5] with the surface bands marked. Here the 
surface bands are identified by applying criterion 1 to the 
1 QL and 2 QL slabs and criterion 2 to the other slabs. 
Critical percentages are self-consistently determined such 
that the wave function projections onto the surfaces are 
stable. The specific critical percentages are as follows: 
45% and 30% for 1 QL and 2 QLs, respectively, and 
60%, 50%, 40%, and 40% for 3 QLs, 4 QLs, 5 QLs, and 
6 QLs, respectively. For very thin slabs (1-2 QLs), cri- 
terion 2 is not applied. When criterion 1 is applied to 
the slabs thicker than 2 QLs, the critical percentages are 
25%, 20%, 16%, and 16% for 3 QLs, 4 QLs, 5 QLs, and 
6 QLs, respectively. Among the surface states marked in 
Fig. [21 only those within 0.2 eV above and 0.3 eV below 
Ef near F survive as surface states. The surface bands 
outside this energy window, far away from F, are strongly 
coupled to the bulk bands [Fig. [Ub)]. A comparison of 
the surface bands identified using the criteria 1 and 2 re- 
veals that the surface states within the energy window 
remain essentially unchanged using the different criteria. 

Now we present a physical interpretation of the surface 
states identified above. Let us first consider a very thick 
slab where the top surface bands do not interact with the 
bottom surface bands. When the thick slab has symmet- 
ric surfaces (the top and bottom surfaces in an identical 
environment), the surface bands become fourfold degen- 
erate at F, because of Kramers degeneracy at the time- 
reversal invariant momentum for both the top and bot- 
tom surfaces. Away from F, the fourfold degeneracy is 
lifted to a twofold degeneracy due to SOC. Neglecting the 
fact that the spin-up and spin-down states are not eigen- 
states in the presence of SOC, for a symmetric thick slab, 
the top surface bands with spin-up (spin-down) are de- 
generate with the bottom surface bands with spin-down 
(spin- up). 

As the slab thickness decreases, the top surface bands 



FIG. 2: (Color online) Band structures of Bi2Te3(lll) films with six thicknesses (surface states marked by the squares). Inset 
in (a): Two-dimensional Brillouin zone. Inset in (b): Schematic view of one QL of a Bi2Te3(lll) film. A, B, and C indicate 
the planar sites of the Te and Bi atoms. 21 = 1.737 A, Z2=2.033 A, 23=2.625 A l2g]. Inset in (c): Zoom-in of the 3 QL band 
structure near Ef = 0. Inset in (e): Zoom- in of the 5 QL band structure near Ef = 0. Inset in (f): A (eV) vs thickness (nm). 



TABLE II: Wave function projections onto each QL (from the 
topmost to the bottommost QL) for the 5 QL slab, as a func- 
tion of energy band and fcy. The specific fcy points for the 
V-shaped (U-shaped) band shown in the inset of Fig. [2je) 
are labeled as Ty, Dv, Ey (Tu, Du, Ey). Each band 
has a double degeneracy. Two wave functions correspond- 
ing to Fv have the following projections onto each QL with 
the opposite spin moments: (0.503,0.204,0.024,0.085,0.189), 
(0.189,0.085,0.024,0.204,0.503). The surface states (SS) are 
identified based on the criterion 2. 





Tu 


Tv 


Dc7 


Gv 


Eu 


Ev 


fc|| (27r/a) 








0.022 


0.022 


0.057 


0.057 


Energy(eV) 


-0.091 


-0.082 


-0.074 


-0.012 


-0.001 


0.111 


topmost 


0.438 


0.503 


0.490 


0.795 


0.328 


0.573 


top— 1 


0.190 


0.204 


0.356 


0.182 


0.323 


0.274 


middle 


0.060 


0.024 


0.128 


0.016 


0.187 


0.101 


bottom-l-1 


0.094 


0.085 


0.023 


0.000 


0.083 


0.038 


bottommost 


0.217 


0.189 


0.003 


0.000 


0.080 


0.017 


SS? 


Yes 


Yes 


Yes 


Yes 


No 


Yes 



interact with the bottom surface bands, resulting in an 
opening of a band gap A, even at T. At F, a linear even 
combination of the top spin-up and the bottom spin-up 



surface states are degenerate with a linear even combina- 
tion of the top spin-down and the bottom spin-down sur- 
face states. At F, the energies of these even combinations 
are separated from those of the similar odd combinations 
(with a double degeneracy) by A. This gap increases 
with increasing fc|| . As the thickness decreases, A in- 
creases. Even though the symmetric slab still possesses 
time-reversal and inversion symmetry, a strong interac- 
tion between the top and bottom surfaces mixes spin-up 
with spin-down states, resulting in a reduction in spin 
polarization at each surface. Thus, spin- flip scattering at 
each surface will no longer be suppressed. However, the 
thin slab must show a remnant of a bulk TI: the number 
of Ef crossings of the surface bands for a given surface 
(either the top or bottom surface). 

Our calculated values of A at F and of the indirect 
band gap Aii^ind corroborate robust surface states for 
three or more QLs. The 1 QL and 2 QL slabs show A 
of 0.4338 and 0.1319 eV at F, respectively. This gap de- 
creases exponentially with increasing thickness and it sat- 
urates at a thickness of 4 QLs [Table HI inset in Fig.[2{f)], 
in agreement with Refs. [19|, |21[. For the slabs four or 
less QLs thick, A does not originate entirely from the 
surface states. At F, for the slabs three or less QLs 
thick, the valence band does not have a surface char- 
acter [Fig. [ira),(b), the V-shaped band in the inset of 



(c)], while for the 4 QL slab, the band slightly below 
the V-shaped band loses its surface character. However, 
as A:|| increases, for the 3 QL (4 QL) slab, the V-shaped 
band (the band slightly below the V-shaped band) re- 
trieves its surface character. The 1 QL and 2 QL slabs 
show AEind of 0.301 and 0.057 eV, respectively. As the 
thickness increases, A_Eind vanishes. 

We now use our accurate identification of surface states 
in the slabs to examine their topological properties. 
For a three-dimensional TI, vq is determined from the 
bulk band structure by a product of parity eigenval- 
ues of all occupied bands (counting Kramers degenerate 
pairs only once) at the eight time-reversal invariant mo- 
menta [l|. However, this procedure is not well defined 
in slabs. Thus, instead, we use an equivalent criterion 
[l|. The surface bands of a TI cross Ef an odd num- 
ber of times between time-reversal invariant momenta 
{vq = 1). With the notations in Ref.[l|], vq can be ob- 
tained from (—1)''" = ''^r^'M^'M '"'m ' where T, M, M2, 
and M3 are four time-reversal invariant momenta in the 
two-dimensional Brillouin zone [inset in Fig. [21[a)]. If the 
surface bands cross Ef an odd (even) number of times 
between F and M, then 7^p7r;|^J = —1 (-1-1). Another 
equivalent criterion states that the conduction and va- 
lence bands with opposite parities are inverted by SOC 
at the Dirac point (F) for slabs with topologically non- 
trivial surface states. These two criteria must provide 
the same result. 

For the 1 QL and 2 QL slabs, the surface states nei- 
ther cross Ef between F and M [Fig. [2ja),(b)] nor be- 
tween M2 and M^. Thus, the 1 QL and 2 QL slabs have 
topologically trivial [vq = 0) surface states. This is also 
corroborated by our finding that the conduction and va- 
lence bands are not inverted at F by SOC. Note that the 
topologically trivial nature of the surface states would 
not change by shifting Ef upward. 

Let us discuss the slabs three or more QLs thick. Their 
surface states are marked in Fig. [2jc)-(f). Within the 
small energy window from Ef, only two surface bands 
are identified. The features of these two surface bands 
discussed below do not depend on the thickness. To bet- 
ter examine the surface states, for example, for the 5 
QL slab, the projections of the wave function onto each 
QL at three fc|| -points for the two energy bands [shown 
in the inset of Fig. [SJe)] a£e hsted in Table HH The V- 
shaped valence band near F is a surface band [Fig. EJc)- 
(f )] . This surface band intersects Ef only once between F 
and M [insets in Fig.[2jc),(e)] and joins the bulk bands as 
fc|| moves further away from F. The occupied U-shaped 
band slightly below the V-shaped band near F [insets in 
Fig. [2Ic),(e)], is again a surface band only within a very 
narrow window in kn , yet it loses its surface character 
before it crosses Ef [see the projections of this band at 
the E point, E[/, in Table [H] and the inset of Fig. Efe)]. 
Additionally, none of the surface bands cross Ef between 



M2 and M3. SOC indeed inverts the conduction and va- 
lence bands. Thus, slabs three or more QLs thick have 
topologically non-trivial (vq — I) surface states, indepen- 
dent of the thickness. Our findings clearly reveal that 
the topological nature of the surface states persists with 
the thickness, differing from previous calculations 21j . 
yet in agreement with experiment [l9|. 

We have constructed a method to accurately iden- 
tify topological surface states within DFT. Using this 
method, we have investigated the topological nature of 
the surface states in thin films of Bi2Te3(lll) using DFT. 
We have found that the topological nature of the surface 
states remains robust with the film thickness and that 
the surface states are topologically protected for films of 
three or more QLs. The method and our findings are 
applicable to thin films of other types of TI. 
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